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Stephanie Jacquot* Benedek Valko^ 



Abstract 

We consider the /3-Laguerre ensemble, a family of distributions generalizing the 
joint eigenvalue distribution of the Wishart random matrices. We show that the bulk 
scaling limit of these ensembles exists for all /3 > for a general family of parameters 
and it is the same as the bulk scaling limit of the corresponding /3-Hermite ensemble. 

1 Introduction 

The Wishart-ensemble is one of the first studied random matrix models, introduced by 
Wishart in 1928 [12]. It describes the joint eigenvalue distribution of the n x n random 
symmetric matrix M = AA* where A is an n x (m + 1) matrix with i.i.d. standard normal 
entries. One can also define versions with i.i.d. complex or real quaternion standard normal 
random variables. Since we are only interested in the eigenvalues, we can assume m + 1 > n. 
Then the joint eigenvalue density on M" exists and it is given by the following formula: 

n J 

r(m-n)-l _lx 



'^n,m+l j<fc k=l 

where /3 = 1, 2 and 4 correspond to the real, complex and quaternion cases respectively and 
Z^m+i is an explicitly computable constant. 

The density ([1]) defines a distribution on for any /3 > 0, n G N and m > n. The 
resulting family of distributions is called the /3-Laguerre ensemble. Note that we intentionally 
shifted the parameter m by one as this will result in slightly cleaner expressions later on. 
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Another important family of distributions in random matrix theory is the Hermite (or 
Gaussian) /3-ensemble. It is described by the density function 

l<i<fe<n fc=l 

on M". For /3 = 1, 2 and 4 this gives the joint eigenvalue density of the Gaussian orthogonal, 
unitary and symplectic ensembles. It is known that if we rescale the ensemble by ^Jn then the 
empirical spectral density converges to the Wigner semicircle distribution ^ — a^^l[-2,2]- 
In [13] the authors derive the bulk scaling limit of the /3-Hermite ensemble, i.e. the point 
process limit of the spectrum it is scaled around a sequence of points away from the edges. 

Theorem 1 (Valko and Virag [I3])- // satisfies n}l^{2^n — |/i„|) — ?■ oo as n —> oo and 

is a sequence of random vectors with density ^ then 

y/in- fil{A^ - /in) ^ Sine/3 (3) 
where Sine^ is a discrete point process with density (27t)^^. 

Note that the condition on yU„ means that we are in the bulk of the spectrum, not too 
close to the edge. The limiting point process Sine/3 can be described as a functional of the 
Brownian motion in the hyperbolic plane or equivalently via a system of stochastic differential 
equations (see Subsection 12.41 for details). 

The main result of the present paper provides the point process limit of the Laguerre 
ensemble in the bulk. In order to understand the order of the scaling parameters, we first 
recall the classical results about the limit of the empirical spectral measure for the Wishart 
matrices. If m/n — 7 G [l,oo) then with probability one the scaled empirical spectral 
measures Vn = }i X]fc=i ^>^k/n converge weakly to the Marchenko-Pastur distribution which is 
a deterministic measure with density 



a\x)= ^^'' f^ i^^lt..,,.](x), a = a(7) = - 1, h = h{^) = I + (4) 

This can be proved by the moment method or using Stieltjes-transform. (See [7] for the 
original proof and [S] for the general /3 case). 
Now we are ready to state our main theorem: 

Theorem 2 (Bulk limit of the Laguerre ensemble). Fix /3 > 0, assume that m/n ^ & 
[1, 00) and let c G (a^, h"^). Let denote the point process given by ([ip. Then 

27ia^{c) (A^ - en) Sine,? (5) 

where Sineg is the bulk scaling limit of the (5-Hermite ensemble. 
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We will actually prove a more general version of this theorem: we will also allow the 
cases when m/n — oo or when the center of the scaling gets close to the spectral edge. See 
Theorem [8] in Subsection 12.21 for the details. 

Although this statement has been known for the classical cases (/3 = 1, 2 and 4) [8], this 
is the first proof for general /3. Our approach relies on the tridiagonal matrix representa- 
tion of the Laguerre ensemble introduced by Dumitriu and Edelman [T] and the techniques 
introduced in [13] . 

There are various other ways one can generalize the classical Wishart ensembles. One 
possibility is that instead of normal random variables one uses more general distributions in 
the construction described at the beginning of this section. The recent papers of Tao and 
Vu [12] and Erdos et al. [3] provide the bulk scaling limit in these cases. 

Our theorem completes the picture about the point process scaling limits of the Laguerre 
ensemble. The scaling limit at the soft edge has been proved in [9], where the edge limit of 
the Hermite ensemble was also treated. 

Theorem 3 (Ramirez, Rider and Virag [9]). If m > n ^ oo then 



where Airy^ is a discrete simple point process given by the eigenvalues of the stochastic Airy 
operator 



Here b'^ is white noise and the eigenvalue problem is set up on the positive half line with 
initial conditions /(O) = 0, /'(O) = 1. 

A similar limit holds at the lower edge: if liminf m/n > 1 then 



Remark 4. The lower edge result is not stated explicitly in [9], but it follows by a straight- 
forward modification of the proof of the upper edge statement. Note that the condition 
liminf m/ra > 1 is not optimal, the statement is expected to hold with m — — )■ oo. This 
has been known for the classical cases /3 = 1, 2, 4 [8]. 

If m — n — > a G (0, oo) then the lower edge of the spectrum is pushed to and it becomes 
a 'hard' edge. The scaling limit in this case was proved in |10] . 
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Theorem 5 (Ramirez and Rider [lO]). If m — n a & (0, oo) then 

where 9^^a is a simple point process that can be described as the sequence of eigenvalues of 
a certain random operator. 

In the next section we discuss the tridiagonal representation of the Laguerre ensemble, 
recaU how to count eigenvalues of a tridiagonal matrix and state a more general version of 
our theorem. Section [3] will contain the outline of the proof while the rest of the paper deals 
with the details of the proof. 

2 Preparatory steps 

2.1 Tridiagonal representation 

In pp Dumitriu and Edelman proved that the /3-Laguerre ensemble can be represented as joint 
eigenvalue distributions for certain random tridiagonal matrices. Let An^m be the following 
n X n bidiagonal matrix: 

X/3(m-l) 

X/3(n-l) X/3(m-2) 

X/3-2 X/3(m-n+l) 

X/3 Xl3{m.-n) _ 

where Xpa^Xph are independent chi-distributed random variables with the appropriate pa- 
rameters (l<a<r?, — l,m- — 1 < b < m — n). Then the eigenvalues of the tridiagonal matrix 
An^mA'^^ are distributed according to the density ([I]). 

If we want to find the bulk scaling limit of the eigenvalues of An^m^^m then it is sufficient 
to understand the scaling limit of the singular values of A^ m-The following simple lemma 
will be a useful tool for this. 

Lemma 6. Suppose that B is an nx n bidiagonal matrix with ai, a2, . . . , a„, in the diagonal 
and 6i, 62, • • • , &n-i below the diagonal. Consider the 2n x 2n symmetric tridiagonal matrix 
M which has zeros in the main diagonal and ai, 61, a2, 62, • • • , On in the off-diagonal. If the 
singular values of B are Ai, A2, . . . , A„ then the eigenvalues of M are ±Aj, i = 1 . . .n. 

We learned about this trick from [2j , we reproduce the simple proof for the sake of the reader. 



Ar 
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Proof. Consider the matrix B 



If Au = \iV and Al^v = XiU then [u, ±v] is 



fi^ 
B 

an eigenvector of B with eigenvalue ±Aj. Let C be the permutation matrix corresponding 
to (2, 4, ... , 2n, 1, 3, . . . , 2r;, — 1). Then C^BC is exactly the tridiagonal matrix described in 
the lemma and its eigenvalues are exactly ±Aj, i = 1 . . .n. □ 

Because of the previous lemma it is enough to study the eigenvalues of the {2n) x [2n) 
tridiagonal matrix 

X/3(m-l) 
X/3(m-l) X/3(n-l) 

Xl3(n-1) X/3(m-2) 



A., 



1 

7^ 



Xl3 





"X.I3{m—n) 





(6) 



X/3(m-n+l) 
X/3 

X/9(m— n) 

The main advantage of this representation, as opposed to studying the tridiagonal matrix 
y4„ m^nm' ^^^^ ^^^^ ^he entries are independent modulo symmetry. 

is an eigenvector for An^m with eigen- 



Remark 7. Assume that [ui, wi, U2, ^2, • • • , w^, u„j 
value A. Then [ui, U2, . . . , u„]"^ is and eigenvector for A'^^An^m with eigenvalue A^ and 



^'1,^2, 



is an eigenvector for A^ m^nm with eigenvalue A^ 



2.2 Bulk limit of the singular values 

We can compute the asymptotic spectral density of An^m from the Marchenko-Pastur distri- 
bution. If m/n — J- 7 G [1, oo) then the asymptotic density (when scaled with ^/n) is 



2\x\a 



J ix - a)ix + a) (b - x)(b + x) ^ .. .. 

— MaMim)- 



(7) 



This means that the spectrum of An,m in is asymptotically concentrated on the interval 
[^/m — y/n, y/rn + y/n\. We will scale around /i^ G {y/m — i/n, y/rn + y/n) where /i„ is chosen 
in a way that it is not too close to the edges. Near /i„ the asymptotic eigenvalue density 
should be close to which explains the choice of the scaling parameters in the 

following theorem. 
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Theorem 8. Fix f3 > and suppose thatm = m{n) > n. Let Kn denote the set of eigenvalues 
of An^m o,nd set 

no = ^na-/"(/i„n-^/2)'-^, n, = n - ^na-/\f^^n-'/')\ (8) 
Assume that as n ^ oo we have 

n/-^nQ^ (9) 

and 

liminfm/ra > 1 or lim m/n = 1 and lim inf i/n > 0. (10) 

n—^oo n— >oo 

Then 

4^/n^{An - fin) ^ Sme^j . (11) 

The extra 1/2 in the definition of no is introduced to make some of the forthcoming 
formulas nicer. We also note that the following identities hold: 

1 2(m + n)^l - (m - nf - /x^ {m - n - filf 
+ 2 4^15 • 4^1 ■ 

Note that we did not assume that m/n converges to a constant or that /i„ = c^/n. By 
the discussions at the beginning of this section (A„, fl ]R+)^ is distributed according to the 
Laguerre ensemble. If we assume that m/n — 7 and /i„ = y/cy/n with c G (0(7)^,6(7)^) 
then both (Q and ffTOj) are satisfied. Since in this case non~^ — )■ d''^{c) the result of Theorem 
IH] implies Theorem |2l 

Remark 9. We want prove that the weak limit of 4y^(A„ — /i„) is Sine/?, thus it is sufficient 
to prove that for any subsequence of n there is a further subsequence so that the limit in 
distribution holds. Because of this by taking an appropriate subsequence we may assume 
that 

m/n — 7 G [1, 00], and if m/n — > 1 then /i„/-\/n c G (0, 2]. (13) 

These assumptions imply that for mi = m — n + ni we have 

liminf mi/n > 0. (14) 

One only needs to check this in the m/n ^ 1 case, when from (fT3!) and the definition of ni 
we get ni/n — )• c > 0. 
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Remark 10. The conditions of Theorem [8] are optimal if hminf m/n > 1 and the theorem 
provides a complete description of the possible point process scaling limits of A^. To see this 
first note that using = (A„ fl M"*")^ we can translate the edge scaling limit of Theorem [3] 
to get 

;(A„- (v^± v^)) ^ ±Airy^. (15) 



If liminf m/n > 1 then by the previous remark we may assume Mm m/n = 7 G (1, oo]. Then 
the previous statement can be transformed into n^^^{An — {\/rn ± ^Jn)) =^ S where H is a 
a linear transformation of Airy^. From this it is easy to check that if n^^n^^ -> c G (0, 00 
then we need to scale A„ — /i„ with "n}/^ to get a meaningful limit (and the limit is a linear 
transformation of Airy^) and if n\^^nQ^ — then we get the bulk case. 

If m/n — 7- 1 then the condition fllU]) is suboptimal, this is partly due to the fact that the 
lower soft edge limit in this case is not available. Here the statement should be true with 
the following condition instead of ( !T0|) : 

2 



IJin\/n{m — n) 



(m — n)^/^ — 7- 00. 



2.3 Counting eigenvalues of tridiagonal matrices 

Assume that the tridiagonal k x k matrix M has positive off-diagonal entries. 

ci a2 62 



M 



Ck-2 



bi > 0, Q > 0. 



If M = [ui, . . . ,Uk] is an eigenvector corresponding to A then we have 



l,...k 



(16) 



where we can we set uq = u^+i = (with cq, hk defined arbitrarily). This gives a single term 
recursion on M U {00} for the ratios 



ro = cx), 



1 



+ A — 



l,...fc. 



(17) 



This recursion can be solved for any parameter A, and A is an eigenvalue if and only if 
Tk = rk\ = 0. 
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We can turn ^ into an angle (pi^x with A — i- (pi^x being continuous, monotone increasing 
and 2 tan(0) = r. Then the values ipkM ^-nd (j)k,Xi will determine the number of eigenvalues 
of M in [Ao, Ai]: 

# {(0fc,Ao,0fc,Ai] n 27rZ} = ^{eigenvalues in (Aq, Ai]} 

This is basically a discrete version of the Sturm-Liouville oscillation theory. 

We do not need to fully solve the recursion f lT7|) in order to count eigenvalues. If we 
consider the reversed version of f[T71) started from index k with initial condition 0: 

rf^ = 0, rt^ = -c,{a,-\ + hrfy\ i=l,...k. (18) 

then A is an eigenvalue if and only if r^^x = ^k-ex- Moreover, if we turn rf^^ into an angle 
0®^ (similarly as before for r and 0) we can also count eigenvalues in the interval [Aq, Ai] by 
the formula 

# {{<PiM - 0f!Ao' </'AAi - 0?aJ n 27rZ} = #{eigenvalues in (Aq, Ai]} (19) 

In our case, by analyzing the scaling limit of a certain version of the phase function (pe^x 
we can identify the limiting point process. This method was used in [13] for the bulk scaling 
limit of the (3 Hermite ensemble. An equivalent approach (via transfer matrices) was used in 
[Ej and [H] to analyze the asymptotic behavior of the spectrum for certain discrete random 
Schrodinger operators. 

2.4 The Sine/3 process 

The distribution of the point process Sine/? from Theorem [1] was described in [13] as a func- 
tional of the Brownian motion in the hyperbolic plane (the Brownian carousel) or equivalently 
via a system of stochastic differential equations. We review the latter description here. Let 
Z he a complex Brownian motion with i.i.d. standard real and imaginary parts. Consider 
the strong solution of the following one parameter system of stochastic differential equations 
for t G [0, 1), A G M: 

dax = :r^=dt + -y^=^ [(e-"^ - l)dZ] , ax{0) = 0. (20) 

It was proved in that for any given A the limit A^(A) = ^ limf^i ax(t) exists, it is integer 
valued a.s. and A^(A) has the same distribution as the counting function of the point process 
Sine/3 evaluated at A. Moreover, this is true for the joint distribution of {N{Xi),i = 1, . . . ,d) 
for any fixed vector {Xi,i = 1, . . . ,d). Recall that the counting function at A > gives the 
number of points in the interval (0, A], and negative the number of points in (A, 0] for A < 0. 
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3 The main steps of the proof of Theorem 



E 



The proof will be similar to one given for Theorem [T] in [13]. The basic idea is simple 
to explain: we will define a version of the phase function and the target phase function 
for the rescaled eigenvalue equation and consider (JTOll with a certain i = £{n). We will 
then show that the length of the interval in the left hand side of the equation converges 
to 27r(A^(Ai) — A^(Ao)) while the left endpoint of that interval becomes uniform modulo 27i. 
This shows that the scaling limit of the eigenvalue process is given by Sine/3. 

The actual proof will require several steps. In order to limit the size of this paper and 
not to make it overly technical, we will recycle some parts of the proof in [TB] . Our aim is to 
give full details whenever there is a major difference between the two proofs and to provide 
an outline of the proof if one can adapt parts of [13] easily. 

Proof of Theorem O Recall that A,„ denotes the multi-set of eigenvalues for the matrix An^m 
which is defined in ([6]). We denote by Nn{X) the counting function of the scaled random 
multi-sets 4nl^'^{An — fin), we will prove that for any (Ai, ■ ■ ■ , A^) G Mf^ we have 

(iV„(Ai), ■ ■ ■ , iV„(A,)) ^ (iV(Ai), ■ ■ ■ , N{\,)) . (21) 

where A^(A) = 2^1imt_j.i ax{t) as defined using the SDE fl20!) . 

We will use the ideas described in Subsection 12.31 to analyze the eigenvalue equation 
A.n,mX = Ax, where x G M^". Following the scaling given in ( ITTl) we set 

A 

In Section H] we will define the phase function Lpi^x and the target phase function ipf-^^ for 
i G [0,no). These will be independent of each other for a fixed i (as functions in A) and 
satisfy the following identity for A < A': 

# {{ipe,x - ^f,x. - ^f,x'] n 27rZ} = iV„(A') - iV„(A). (22) 

The phase function ip will be a regularized version of the phase function obtained from the 
ratio of the consecutive elements of the eigenvector. The regularization is needed in order 
to have a phase function which is asymptotically continuous. Indeed, in Proposition [16] of 
Section [5] we will show that for any < £ < 1 the rescaled version of the phase function 
(fie^x in [0,?^o(l — converges to a one-parameter family of stochastic differential equations. 
Moreover we will prove that in the same region the relative phase function ag^x = fe,\ — y^e,o 
will converge to the solution ax of the SDE (!20|) 



aLno(i-£)J,A =^ ax{l -e), as n 00 (23) 
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in the sense of finite dimensional distributions in A. This will be the content of Corollary [T71 
Next we will describe the asymptotic behavior of the phase functions yj^^A, «£,a and y^®^ 
in the stretch i G [[^0(1 — e)\, 77-2] where 

n2=[no-}C{n\/^yi)\. (24) 

(The constants e, K, will be determined later.) We will show that if the relative phase function 
is already close to an integer multiple of 27r at [r2o(l ^ ^)J then it will not change too much 
in the interval [[?t-o(1 — £^)J;'^2]- To be more precise, in Proposition [THl of Section [B] we will 
prove that there exists a constant c = c(A, /3) so that we have 

E [\{a\no{i-e)\,\ -ttna.A) A l] < c dist(aL„o(i-£)j,A, 27rZ) + v/e + no^^^(ni^^Vlogno) + K,^^ 

for all /C > 0,e G (0,1),A < |A|. 

We will also show that if /C — t- 00 and /C(n|^^ V V)nQ^ — )■ then the random angle v^n2,o 
becomes uniformly distributed module 27r as ri — )■ 00 (see Proposition l22l) . 

Next we will prove that the target phase function will loose its dependence on A: for 
every A G M and /C > we have 

«®A = V%,x - fn2,o -^0' as n ^ 00. (25) 

This will be the content of Proposition [23] in Section [71 

The proof now can be finished exactly the same way as in [13]. We can choose e = 
e{n) and K, = lC{n) — )■ 00 so that the following limits all hold simultaneously: 

{a\no(i-e)\,\,i = ^ ■ ■ -d) ^ {2TiN{\i),i = l...d), 

^n2,o — ^ Uniform[0, 2ti] modulo 27r, 
P 

«L"o(l-£)J,Ai - «n2,A, > 0, 2 = 1,...,^, 

«n2,A, 0, i = l,...,d. 

This means that if we apply the identity (l22l) with A = 0, A' = Aj and i = n2 then the length 
of the random intervals 

converge to 27riV(Aj) in distribution (jointly), while the common left endpoint of these inter- 
vals becomes uniform modulo 2n. (Since v'n2,o and ip^^ q are independent and (pn2,o converges 
to a uniform distribution mod 27r.) This means that 7^{2/c7r G Jj : A; G Z} converges to N{Xi) 
which proves ( 12T|) and Theorem [8] □ 
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4 Phase functions 



In this section we introduce the phase functions used to count the eigenvalues. 



4.1 The eigenvalue equations 



Let Sj = y/n — j — 1/2 and pj = — j — 1/2. Conjugating the matrix An,., 
{2n) X (2n) diagonal matrix D = D^") with diagonal elements 



with a 



D2i,2i = n 



Xl3{m-l)Xl3{n~e) 



D 



2i+l,2't+l 



T ~X.I3{m—t)"X.li{n- 

~7Fp 



n 



we get the tridiagonal matrix A^^ = D ^An^mD: 



A, 



D 



Po + Xo 
Pi So + Yq 



Sl 







Pi + Xi 



Pn^l 



where 



X/3(m-i-l) 



(26) 



Pp. 



Pe, 0<i<n-l, 



Sn-2 + Yn-2 
Sn-1 Pn^l + Xn^l 

Pn 



^+1 



The first couple of moments of these random variables are explicitly computable using the 
moment generating function of the ^^-distribution and we get the following asymptotics: 



= C((m 



-3/2^ 



, EX| = 2//3 + C((m EXf = 0{l), 



EF^ = C((n-£)-3/2)^ EF/ = 2//3 + C((n EF/ = C(l), 



(27) 



where the constants in the error terms only depend on /3. 

We consider the eigenvalue equation for A^^ with a given A G M and denote a nontrivial 
solution of the first 2n — 1 components by ui, vi, U2, V2, ■ ■ ■ , Un, Vn- Then we have 



Pe+iue+i + {se + Ye)ui+2 = Af^+i, 



0<i<n-l, 
0<i<n-2, 
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where we set Vq = and we can assume Mi = 1 by linearity. We set = r^ A = ui+i/vi, 
< i < n — 1 and = f^^A = ^i/ui, 1 < i < n. These are elements of M U {c>o} satisfying 
the recursion 



with initial condition tq = oo. We can set y„ = and define r„ via f l2^ with £ = n — 1, then 
A is an eigenvalue if and only if r„ = 0. 

4.2 The hyperbolic point of view 

We use the hyperbolic geometric approach of [13] to study the evolution of r and f. We will 
view MU{oo} as the boundary of the hyperbolic plane EI = {'^z > : 2 G C} in the Poincare 
half-plane model. We denote the group of linear fractional transformations preserving EI by 
PSL(2,M). The recursions for both r and f evolve by elements of this group of the form 
X ^ b — a/x with a > 0. 

The Poincare half-plane model is equivalent to the Poincare disk model U = {\z\ < 1} 
via the conformal bijection lJ{z) = ^ which is also a bijection between the boundaries 

= M U {oo} and 9U = {|^| = 1, ^ G C }. Thus elements of PSL(2, M) also act naturally 
on the unit circle dU. By lifting these maps to M, the universal cover of dU, each element 
T in PSL(2,M) becomes an M M function. The lifted versions are uniquely determined 
up to shifts by 27r and will also form a group which we denote by UPSL(2,M). For any 
T G UPSL(2,R) we can look at T as a function acting on dM, dV or R. We will denote 
these actions by: 

dM^dM : z.T, dV ^ dV : z ^ ZoT, ^ : z ^ 2*T. 

For every T G UPSL(2, M) the function x t— /(x) = x*T is monotone, analytic and quasiperi- 
odic modulo 27r: /(x + 27r) = /(x) + 27r. It is clear from the definitions that e*^oT = e*-^'-^'-' 
and (2tan(x)).T = 2tan/(x). 

Now we will introduce a couple of simple elements of UPSL(2, M). For a given a G M we 
will denote by Q(a) the rotation by a in U about 0. More precisely, (/5*Q(a) = ip + a. For 
a > 0, 6 G M we denote by A(a, h) the affine map z — ?■ a{z + h) in EI . This is an element of 
PSL(2,R) which fixes oo in H and -1 in 9U. We specify its lifted version in UPSL(2,M) by 
making it fix tt, this will uniquely determines it as a M — M function. 




(28) 



(29) 
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Given T G UPSL(2, M), x, y G M we define the angular shift 

ash(T, X, y) = {y*T — x*T) — {y — x) 

which gives the change in the signed distance of x, y under T. This only depends onv = e*^, 
w = and the effect of T on dU , so we can also view ash(T, ■) as a function on dU x dU 
and the following identity holds: 

ash(T,w,M;) = argrg JtyoT/t^oT) - argrg o^Jui/t;). 



The following lemma appeared as Lemma 16 in [13], it provides a useful estimate for the 
angular shift. 

Lemma 11. Suppose that for a T G UPSL(2,]R) we have {i + z).T = i with \z\ < ^. Then 



ash(T, V, w) 



3f? 



_ i{2+v+w) 2 
^ 4 ^ 



[W — V 

= -"^Kw - v)z] + 62 = Si, 

where for d = 1, 2, 3 and an absolute constant c we have 

kdl < c\w - v\\zf < 2c\zf. 
If V = —1 then the previous bounds hold even in the case \z\ > |. 

4.3 Regularized phase functions 

Because of the scaling in (fTTl) we will set 



^3 



(30) 



(31) 



4nn 



We introduce the following operators 



J, = Q(7r)A(s,/p,, fin/se), M, = A((l + X,/p,)"\ A/(4ny V))A(-^, 0), 



1/2 A / Pi 

J, = Q(7r)A(p,/s,, ^ln/Pi). M, = A((l + Y,/s,Y\ \/{Anl'%,)). 

Then (128|) and ( 129|1 can be rewritten as 

r^+i = r^J^M^^M^, Tq = oo. 

(We suppressed the A dependence in r and the operators M, M.) Lifting these recursions 
from dW to M we get the evolution of the corresponding phase angle which we denote by 

= (f)e*JiMeJeMe, 0o = -tt. (32) 
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Solving the recursion from the other end, with end condition we get the target phase 
function ^^^^ ■ 



bf = 0f+,*M7iJ7^M7^J7\ 0« = 0. (33) 

It is clear that and (pf^ are independent for a fixed £ (as functions in A), they are 
monotone and analytic in A and we can count eigenvalues using the formula (fT9|) . 

Note that and do not depend on A and they are not infinitesimal. The main part 
of the evolution is J^J^. This is a rotation in the hyperbolic plane if it only has one fixed 
point in H. The fixed point equation p£ = p^.J^J^ can be rewritten as 

se{pe - pese) 

This can be solved explicitly, and one gets the following unique solution in the upper half 
plane if i < Uq: 

fil-m + n . I {pl~m + n)^ 

(One also needs to use the identity pj — sj = m — n.) This shows that if £ < no then J^J^ 
is a rotation in the hyperbolic plane. We can move the center of rotation to in U by 
conjugating it with an appropriate affine transformation: 

Jfj£ = Q(-2arg(p^pf))^^ \ 

Here = A($5(p£)^\ -3fJp^), = Y-^XY and 



pl + m-n . / (p^ + m - 
In order to regularize the evolution of the phase function we introduce 

where = Q(2 arg(poPo)) • • • Q(2 arg(p£p^)) and Q_i is the identity. It is easy to check 
that the initial condition remains (f^^x = vr. Then 

= ife*Qi\Tj'Q{-2 arg(p,))^"M^^M,T,T7iT,+iQ, 



Note that the evolution operator is now infinitesimal: M^, and ^T^+i are all asymp- 
totically small, and the various conjugations will not change this. 
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We can also introduce the corresponding target phase function 



^tx ■= 0?A*T^Q^-i. < £ < no. (36) 

The new, regularized phase functions ip(,^\ and ip^y^ have the same properties as (f),4>^, i.e.: 
they are independent for a fixed i (as functions in A), they are monotone and analytic in A 
and we can count eigenvalues using the formula 

We will further simplify the evolution using the following identities: 



— + b 
r \ a 

From this we get 
where 



r — 6 Q ar 



b-i 
b + i 



r.JeTi 



I Pe ^ /in 



J,T, = t,Q,(-2arg(p^)) 



This allows us to write 



- 2- — pn 



^pm 



^Pi 



r 1 

^PePeSeJ se^pi ^pe 



'^Pi '^pi 



Mi 



(m; 



TMQ(-2arg(p,))Qf 



(m; 



(37) 



where 



Thus 



= QfQ(-2arg(p^)). 



Qt 



We will introduce the following operators to break up the evolution into smaller pieces: 



L,,, = A(l, X/iAnl^'p,)), L,,A = A(l, X/i4n'/'s,)) 
Se,x = Lji (tj'A{^{1 + X,/pe)-\ 0) t, 
S,,A = LS (Tr'A((l + Ye/s,r\ 0) T,+i) . 



1/2, 



Then 



ift+i = (pi*{L/] (S^^o) L/ 8^,0 = V^^* (S^a) S^,A 
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5 SDE limit for the phase function 



Let Ti denote the a-field generated by ^j,x,j < i- Then (pi^\ is a Markov chain in i 
with respect to We will show that this Markov chain converges to a diffusion limit 
after proper normalization. In order to do this we will estimate K[lp£^i^x — Lpi^x\J^e] and 
E [{ipi+i^x — f£,x){f£+i,x' — 'fi,x')\J^£] using the angular shift lemma, Lemma [TTl 



5.1 Single step estimates 

Throughout the rest of the proof we will use the notation k = no — i. We will need to 
rescale the discrete time uq in order to get a limit, we will use t = i/no and also introduce 
s{t) = ^/l — t. We start with the identity 



V 2 



= \/nQ- i = \fk= ^JrTosit) . 
Note that this means that 



n-no-1/2 no-i / rii k 



. _ j m - no - 1/2 ^ . / riQ- i _ j mi j k 

~ \l m-£-l/2 m-i-l/2~ y mi + k^l mi + k ^^^^ 



where the sign in 3ftp£ is positive if pn > y/m — n and negative otherwise. 
For the angular shift estimates we need to consider 



Zi,x = ^-S^^;, - t = , . + -- — r-T + ^7 ] =Vi + ve, 

^s{t) pe V 4nos(t) Pe^pe 



v • c-i • P^^^ < f -^1 Pi+i ~ Pi \ T> I - un\ 



We have the following estimates for the deterministic part (by Taylor expansion): 



where = p^'^'{t) = ^Jm/nQ — t and p(t) = p^'^'{t), p{t) = ff'{t) are defined by equations 
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(!38|) and ( 139|) with £ = not. For the random terms from (127|) we get 
El^^ = Ef^ = C(A;-i/2^n-£)-3/2), 

E^/ = Lq{^\t) + 0{k-\n-i)-^), EV,^ = —q^^\t) + 0{k-\n 
no no 

E\V,^\ = E\V,^\ = —q^''\t) + Oik-\n-£y^), E\V.^\,E\%^\ = Oik-'^^^), d = 3,4, 
no 

where the constants in the error term only depend on /3 and 
We introduce the notations 

and Ai/2fx,\ = fx+y2,x - fx,x, Afx,x = fx+i,x - fx,x- We also set for i e Z+ 

2 -s2 2 -^2 2-^2 

Ve = PoPoPiPi ■ --PiPi- 

Remark 12. We would like to note that the 'half-step' evolution rules ipe^\ — )■ (/9£+i/2,a, 
</5£+V2,A — ^ '^e+i,x are very similar to the one-step evolution of the phase function ip in |13j . 
Besides the fact that the £ — > £ -|- 1/2 and £ + 1/2 — j-f + l steps are quite different from each 
other the other big difference between our case and [13] is that here the oscillating terms 
Qi, Qi are more complicated. 

The following proposition is the analogue of Proposition 22 in [13] . 

Proposition 13. For i < no we have 

E [Av.mlv^^A = x]= -b^^\t) + -osc(^) + 0{k-'/^) = 0{k~') 
' no no 

E [Ai/2(pe.xAi/2Vi,x'\<^e,x = x,(pe^x' = y] = — a''^\t,x,y) H osc^^^ + 0{k~-^^'^) 

no no 

E [Av,<^,h_v.,a|(^£+V2,a = x\= -hf{t) + -osc(3) + 0{k-^'^) = 0{k~^) 

no no 

E [Ai/2y5£+i/2,AAi/2V?£+i/2,A'|</3f+V2,A = X, v?£+i/2,A' = v] = — a^^\t , X , y) H osc^^^ + 

no no 

E[|Av2¥^£,a|>aa = x] ,E[\A,/,^e+i/2,x\''\fi,x = x] =0{k-'"^), d = 2,3 
where t = i/no, 



, (1) _ A ^ M ^ /,(2) = A _ ^ 

"A 4s^2p^^2^' "A 4^ ^ 2(3s 



2 ■ 
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The oscillatory terms are 



osc 



(1) = g?((-t;,-Vi)/2)e-^^p72^,) + 3?(ze-2-p7%2^«)/4, 



osc 



Proof. We start with the identity 

</'£+V2,A-</2£,A = </2^+i,A*Q7^-v^^,A*Q7^ = '/'^,A*Q7^S^,A-V5£,A*Q7^ = ash(S^,A, e*'^''^w7^ 

Here we used the definition of the angular shift with the fact that S^^a (and any affine 
transformation) will preserve oo e H which corresponds to —1 in U. A similar identity can 
be proved for ^1/2^1+^/2, \- 

The proof now follows exactly the same as in [TB], it is a straightforward application of 
Lemma [TT] using the estimates on v^^\, Vi^x, Vi, Vi. □ 

5.2 The continuum limit 

In this section we will prove that Lp^'"'\t, \) = fitno\,x converges to the solution of a one- 
parameter family of stochastic differential equations on t G [0,1). The main tool is the 
following proposition, proved in [13] (based on [TT] and [1]). 

Proposition 14. Fix T > 0, and for each n > 1 consider a Markov chain X" G M.'^ with 
i = 1 . . . lnT\ . Let ^^"(x) be distributed as the increment — x given = x. We define 

b^{t,x) = nE[yL:,j(x)], a-{t,x) = nE[YC,,^{x)Y;:,,^{xn 

Suppose that as n —)■ 00 we have 

|a'^(t,x)-a"(t,|/)| + |&"(t,x)-6"(t,y)| < c\x - y\ + 0(1) (42) 

supE[\Yp{x)f] < cn^l\ (43) 

and that there are functions a, b from M x [0, T] to M'^'^, M'^ respectively with bounded first and 
second derivatives so that 



sup 



a^{s, x) ds 













/ a{s,x)ds 


+ sup 


[ h'\s,x)ds 


Jo 




Jo 



Assume also that the initial conditions converge weakly, Xq =^ Xq. 



Then (-^^^j , < t < T) converges in law to the unique solution of the SDE 

dX = bdt + (TdB, X(0)=Xo, t e [0,T], 

where B is a d-dimensional standard Brownian motion and a -.W^ x [0, T] is a square root 
of the matrix valued function a, i.e. a(t,x) = a{t,x) a(t,x)'^ . 

We will apply this proposition to ipe^x with i < no{l — e) and £ G Z/2, so the single steps 
of the proposition correspond to half steps in our setup. 

The following lemma shows that the oscillatory terms in the estimates of Proposition 
[T5] average out in the 'long run'. Its proof relies on Proposition [T3] and Lemma [23 of the 
Appendix. 

Lemma 15. Let |A|, |A'| < A and e > 0. Then for any ii < no(l — e), £i G Z 



- V E [Ay, ^e,x I ^e,x = x] = - V ^^(t) + 0{n-'/' + n'/^'^') (45) 

nn Tin ^ — ' 

1 ~ 1 ^'"^ 

— E [Ai/2 V5£,aAi/2 v3£,a' I = X, V2£,A' =y] = — a{t, x, y) + 0{nQ^'^ + n^'^n^^''^) 

^ 0<f<£i ^ £=0 

where t = i/uQ, the functions bx, a are defined as 

anc? i/ie implicit constants in O depend only on e, f3,X. The indices in the summation ^ run 
through half integers. 

Proof of Lemma [75l We will only prove the first statement, the second one being similar. 
Note that bx{t) = b'^^\t) + fef (t). 

Summing the first and third estimates in Proposition [T2] we get fl45p with an error term 

. h-i ii-i 

- V 3fJ(ei,,r^,) + - V 3fi(e2/r/f) + 0{n-'/'), (47) 

where the first two terms will be denoted Ci;C2- Here 

ei,e = {{-vx - iq^'^/2)p1 + {-Vx - iq^'^ /2)) e-'^ 62,. = i{pj'q^''> + q^'^)e-''y4 

where for this proof c denotes varying constants depending on e. Using the fact that 
Vx, Vx, q^^\ q^"^^ and their first derivatives are continuous on [0, 1 — e] we get 

ki/l < c, \ei^e - Ci^e+il < cng \ (48) 
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Applying Lemma [25] of the Appendix to the first sum in (H7|l : 

1 1 ^'"^ 

no no 

Since ii<no{l-e)we have jF^^^ | < c(l + ni^^/^-^/^) < c(ny^no + 1) and 

ICil < ciriQ^^^ny^ + no^). 
(Recall that k = rio — i.) For the estimate of ^2 we first note that 

I I 1 ^0 I --2 I 2| 1 ^0 I -2 2 , 1 I /An\ 

M = ^^|P, +P.I = ^x"^^^^ + '"- 

We will use Lemma [25] if |p|p^ + 1| is 'big', and a direct bound with f H9|) if it is small. To 
be more precise: we divide the sum into three pieces, we cut it at indices and £2 so that 

Ipip'j + 1| < n^^^'^ if A; G [A;2, kl] and |p^p^ + 1| > n^^^'^ otherwise. (50) 

Note that one or two of the resulting partial sums may be empty. We can always find such 
indices because axg pjpj is monotone if /i„ > i/m — n and if /i„ < y/m — n then argpfpf 
decreases ii k > y/m^n^ then it increases. (See the proof of Lemma [25]) 

We denote the three pieces by ^2,4, « = 1, 2, 3 and bound them separately. Since k > eno, 
Lemma [25] gives 

l> I ^ / 1/2 -3/2 , -1/2n 

|C2,i| < c(n/ Uq ' ' ). 

The term 1^2,31 can be bounded exactly the same way, so we only need to deal with (2,2- Here 
we use ( 149|) to get a direct estimate: 



1 1 

k£[k2 ,k^]nleno ,no] 



~2 , 2| ^ -1/2 

+ Pi\ <cno'. 



Collecting all our estimates the statement follows. □ 

Now we have the ingredients to prove the continuum limit. 

Proposition 16. Suppose that m/no — ?■ k G [l,oo], n/rio — )■ z/ G [l,C)o] and that eventually 
fin > \Jni — n or < \/m — n. Then the continuous functions p^"'\t)~^ , p^''^\t) , p^"'\t) 
converge to following limits on [0, !).■ 



p-\t) = {K-t)-^/\ p(t) = ±J—-+zJ-—. p{t) 



u — t ' y u — t' ' y K — t ' y K — t^ 
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where the sign in depends on the (eventual) sign of — \/m — n. If k = oo then 
p~^{t) = 0, p(t) = 1 and if u = oo then pit) = ±1. 

Let B and W he independent real and complex standard Brownian motions, and for each 
A G M consider the strong solution of 



dip 



A 



<^a(0) = TT. (51) 



Then we have 



V\\not\ =^ as n OO 



where the convergence is in the sense of finite dimensional distributions for A and in path- 
space D[0, 1) for t. 

Proof. The proof is very similar to the proof of Theorem 25 in [T3]. One needs to check that 
for any fixed vector (Ai, . . . , A^) the Markov chain (v5^,a,! 1 < < c^), ^ < \ ^)'^oJ , ^ G Z/2 
satisfies the conditions of Proposition [H] and to identify the variance matrix of the hmiting 
diffusion. Note that because our Markov chain hves on the half integers one needs to slightly 
rephrase the proposition, but this is straightforward. 

The Lipshitz condition f H2]) and the moment condition fH51) are easy to check from Propo- 
sition [131 The averaging condition ( 14^ is satisfied because of Lemma [151 using the fact that 
because of the conditions of the proposition, the functions b^(t),a(t,x,y) converge. This 
proves that the rescaled version of [ipe^x^, 1 < J < c?) converges in distribution to an SDE in 
M*^ where the drift term is given by the limit of {b^,j = 1 . . .d) and the diffusion matrix is 
given by a{t,x)j,k = [e^(-'=-^)] + 

The only step left is to verify that the limiting SDE can be rewritten in the form f[5T]) . 
This follows easily using the fact that if Z is a complex Gaussian with i.i.d. standard real 
and imaginary parts and uji,uj2 E C then 

The following corollary describes the scaling limit of the relative phase function ae^x. 

Corollary 17. Let Z be a complex Brownian motion with i.i.d. standard real and imaginary 
parts and consider the strong solution ax{t) of the SDE system i^) . Then tt[notJ,A ==^ o:x{t) 
as n 00 where the convergence is in the sense of finite dimensional distributions for A and 
in path-space D[0, 1) for t. 
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Proof. We just need to show that for any subsequence of n we can choose a further sub- 
sequence so that the convergence holds. By choosing an appropriate subsequence we can 



assume that m/riQ, n/uQ both converge and that /i^ — y/m — n is always positive or nonneg- 
ative. Then the conditions of Proposition [16] are satisfied and ax = ^Px — will satisfy the 
SDE ( !20l) with a complex Brownian motion Zt := 6^'^°^^'' dWt- From this the statement of 
the corollary follows. □ 

6 Middle stretch 

In this section we will study the behavior of ae,x and ipe^x in the interval [[(1 — e)nQ\,n2] 



with 77-2 



no - IC{n/\l) 



. The constant /C will eventually go to oo, so we can assume 



that /C > Co > with Cq large enough. 
6.1 The relative phase function 

The objective of this subsection is to show that the relative phase function ai^x does not 
change much in the middle stretch. 

Proposition 18. There exists a constant c = c(A,/3) so that with y = rig ^''^(n^^^V log^o) 
we have 

E [\{ae,,x - ae,,x) A 1| J-,J < c {d{ae,,x, 27iZ) + + y + /C-^) (52) 
for a// /C > 0,e G (0, 1), A < \X\,no{l-e) < /i < /2 < ^2, £ G Z. 

Because of the moment bounds (1271) we may assume that 

l^^l, \Yi\ < ^V^s{i/no), for £ < n^. (53) 

Indeed, the probability that (135]) does not hold is at most c(no — n2)~^ < c/C^^ which can 
be absorbed in the error term of (1521) . 

We first provide the one-step estimates for the evolution of the relative phase function. 

Proposition 19. There exists c = c{(3, A) so that for every £ < n2 and \X\ < \ we have the 
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following estimates 

riQ 

Uq 

= Oia,,,k-' + k-'/^^') (54) 
E (Aal Jj-,) = 0(d,,AA;-' + /c^'rio 1) (55) 
E (I Aa,,AA^,,A| I J-^) = 0{ai,xk-') (56) 

where ai^\ denotes the distance between ai^\ and 2n. 

Proof. We first prove estimates on Ai/aO^^A and Ai/2a^+i/2,A- In order to do this, we break 
up the evolution of ipi^x into even smaller pieces: 

where 77^+1/4, a and v2£+3/4,a are defined accordingly. We also define the relative phase functions 
for the intermediate steps in the natural way. 

By choosing c(/3, A) large enough we can assume j^^^ < ^ for £ < ^2 < n — /C. Using 
this with the cutoff (153|) the random variables Zi^x, Zg^x defined in ( HOj) are both less than 
1/3 in absolute value. This means that we are allowed to use Lemma [TT] in the general case 
for each operator appearing in fl57p (i.e. the condition \z\ < 1/3 is always satisfied). From 
this point the proof is similar to the proof of Proposition 29 in [13] . We first write 

Ay, ae,x = ash(L7^ -1, e'^^'-'vepf) + a.sh{S%, e'^'+y^'^vepf, e'^'^'^wf) 

+ ash(S,^o^ e'^^'-^wf, e'^'^'^wf). (58) 

Using Lemma [TT] one can show that the first two terms in 0581) are 

of C»(no^/^A;-i/2). Using 

Lemma [TT] again for the third term together with 

jgic^f,^ _ giv'f.oi = \^i^i.x _ 1| < a£,A, |e'^'^^'^ - e'^'^^-^l < 2d£,A 
we get the analogue of f[M|) for Ai/2a£,A: 

E (Av.«^,a|-F^) = --3ft [Vi (e"''"^'^ - e-^^^'°) [pf [vx + iq^'^ /2) ]} 

- -3ft - e-2^^^-°) [/57^g(^) ]} + 0(«.,aA;-^/^ + k-^'W") 

= 0{ae,xk-^ + k-^^^n'^^^) (59) 
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We can prove the analogues of (155|) and (156!) and similar bounds for Ai/2a£+i/2,A the same 
way. We can also prove 



(60) 



this is the analogue of Lemma 32 from [13] and it can be proved exactly the same way. 
To get ( l54l) we write 



E [Aat,\\ipifi = X, ipi^x = y] 

= E [Ai/2a£,A|v5£,o = x,(pi^x = y] + E [E[Ai/2ae+i/2,\\J^£+i/2]\ve,o = x,(pe^x = y] 
= E [Ai/^ae^x\vifi = x,(pi^x = y] +E [Ai/2ai+i/2,x\^e+i/2fl = x,(pi+i/2,x = y] 
+0{&,^xk-^ + k-''/^n-"^) 

where the last line follows from fl60|) and the just proved half step estimates. Now applying 
fl59|) and the corresponding estimate for Ai/2a^+i/2,A "we get fl5^ . The other to estimates 
follow similarly. □ 

The next lemma provides a Gronwall-type estimate for the relative phase function. This 
will be the main ingredient in the proof of Proposition [181 The proof is based on the single 
step estimates of Proposition UHl and the oscillation estimates of Lemma l25l the latter will 
be proved in the Appendix. 

Lemma 20. There exist constants Cq,Ci,C2 depending on X,I3 and a finite set J depending 
on n, Hi, Till so that with y = ^(^^^^(n}^^ V log no) we have 

i2-2 

|E(a^2,A - ae^,x\J^e,) \ < Ci{y + y/e) + E(a£2-i|-^4)/2 + ^ 6fE(d;£| J"^) 

0<be< Co{n/^k^^/^ + A;-3/2 _^ ^^Six{n/\ l)k^'- + l(^6j)) 
if IC> C2, |A| < A and no(l - e) < ii < £2 < ^2- 

Proof. Recall that k = Uq — i, ki = uq — ii. We denote xi = E[d^|J-'^J and set 
no 

e2,e = --(e-2*^^-^-e-2^^^>o)[/)7V'^ + g^'V4. 



From Proposition UHl we can write 



|E[a2-«i|J-,J| < 



e=ei 



+ 



£2-1 



£2-1 



£=£1 
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whose terms we denote Ci; C2, Cs and (4 respectively. Clearly, (3 is of the right form and 

|C4|<£a:-VV/'<cv/F, 

k=l 

SO we only need to bound the first two terms. 
We will use 

which is the 'one-step' version of flBU]) and can be proved the same way as Lemma 32 in 
From this we get the estimates 

\ei,e\ < cxe/k, |Aej_£| < ck~'^xi + criQ^^'^k'^^'^ . 

Then by Lemma [25] we have 

with Fi^^i < C{n\^'^k~^/'^ + 1). Collecting the estimates and using k2 > /C(n^^^Vl) we get 

£2-2 

ICil ^ clC~^X£2-i + c?7,Q ^^^max(n}^^, 1) + Xf {ny^k^^^'^ + /c^^). 

e=ei 

In order to bound (2 we use a similar strategy to the one applied in the proof of Lemma 
We divide the index set [i*!, ^2] into finitely many intervals Ii, I2, ■ ■ ■ , la so that for each 
1 ^ J ^ o one of the following three statements holds: 

for each i G Ij we have k > ^Jn^rn[ and + 1| > k"^^"^, (61) 
for each I G Ij we have k < ^Jri^m{ and \p\p\ + 1| > k^^^"^, (62) 
for each £ G Ij we have \plp1 + 1| < k'^^"^ . (63) 



/2 



It is clear that if we divide [£i, £2] into parts at the yjnimi and the solutions of \pip1 + 1 
A;~^/^ then the resulting partition will satisfy the previous conditions. Since |p^p| + l| = k~^ 
has at most three roots (it is a cubic equation, see the proof of Lemma |25] for details) we 
can always get a suitable partition with at most five intervals. Moreover the endpoints of 
these intervals (apart from £1 and ^2) will be the elements of a set of size at most four with 
elements only depending on n,mi,ni. 
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We will estimate the sums corresponding to the various intervals Ij separately. If Ij 
satisfies condition (jHS]) then we use 

\e2A = ^\e-'- - e-'^y\l\py, + 1| < ck~'/'x, 

to get 

\J2'^{Vie2,i)\<cY,k-'/'x, 

If Ij = satisfies condition (1^ then we use Lemma U5\ to get 

q-i 

e&ij i=ii 
where we have 

l^iS/jl < I^Sl + 1 < c{k'" + nT{kl)-'" + 1). (64) 
We can bound the first term as 

xt^^{kl)-\{klf/^ + n^\k*^r^'^ + 1) < c/C"i/2^,. 

using k > max{ni , 1). For the general term in the sum we get 

c(fci/2 + ny\k;r'/' + l){xek-' + n-"\~''^) 

< cx^(A;-3/2 ^ Jnax(ni/^ l)k-^) + c{n-^'^k-^ + max(nl/^ l)n-^'^k-'^'^) (65) 

Note that the sum of the error terms fl65l) is 

c (no ^^"^k'^ + max(ni^^, l)nQ "'^^^/c"'^/^) < c(nQ ^^"^ log no + max(ny^, l)nQ ^^^) 

where we used m.a,'x.{n^ , 1) < < nQ. Putting our estimates together: 

I 3?(r/|e2,£) I <c/C-i/2^,. + c ^ x,(A;-=^/2 ^ inax(ni/', 1) A;-^) (66) 



+ c(ng log no + max(nj^^, l)nQ 



The only case left is when Ij = [i^, i*^ satisfies condition (1621) . If /i„ > -^/m — n then we have 
the same estimate for F^^! as in ( |64l) so we get exactly the same bound as in ( l66l) . If we 
have /i„ < -y/m — n then we use (17^ of Lemma (l25il with the bound 

|i^^il<c(A;^/^ + nf(A;,rV2 + i). 
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Copying the previous arguments we get 

l&Ij £=£l 

I /■ -1/2 1 , ^ 1/6 1 \ -1/2n 

+c[nQ logrio + niax^ra^' Aj^o )■ 

Collecting our estimates, noting that ^2 ~ 1 is the endpoint of one of the intervals Ij and 
letting /C be large enough we get the statement of the lemma. □ 

The proof of Proposition [18] relies on the single step estimates of Proposition [19] and the 
following Gronwall-type lemma which was proved in [T3] . 

Lemma 21. Suppose that for positive numbers X£, bi, c, integers ii < i < £2 we have 

e-i 

xe < ^ + c + b^Xy (67) 

Then xi^ < 2 {xi^ + c) exp ^3 Ylf=el ^j) ■ 

Now we are ready to prove Proposition [T8l 

Proof of Proposition [Tgi We will adapt the proof of Proposition 28 from [13]. Let a = a^^,A 
and define a(),a^ G 27rZ so that [ao,a^) is an interval of length 27i containing a. We can 
assume that A > 0, the other case being very similar. We will drop the index A from a and 
we will write E(.) = 'K{.\J^i^). 

We will show that there exists cq so that if /C > cq, then if a = or then 

E|a^2 — a| < ci{\a — d\ + y/e + y) . (68) 

The claim of the proposition follows from this by an application of the triangle inequality, 
the additional condition k > cq is treated via the error term 

In order to prove ( [681) for a = we follow the steps described in Proposition 28 from 
|13] . Using the exact same argument we only need to prove that for the coefficients bi in 
Lemma [20] are bounded by a constant depending only on A, (3 and that a never goes below 
an integer multiple of 27r that it passes. 

The first statement is easy to check, we have, 

'f , , . 1, . 1) . ..ax(„:^ 1, !).«.))<.. 

i=ll 
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To prove the other statement first recall the evolution steps ( 157|) and that aj^x = fj^\ — (Pj,x for 
j G Z/4. Using the fact that the maps L^ ;^; ^ and their conjugates are monotone in A (as 

/4,A and a£+i/2,A < "£+3/4 



functions on M) we get that ag^x < tt£+V4,A and ai+ih^x < tt£+3/4,A- Moreover, since (S^^o)'"^ 



and ( S^^o) are 27r— quasiperiodic functions on M we have ,aJ 27r = ['^£+1/2 ,aJ 271 and 

Lq!£+3/4 ,aJ 27r = L'^£+i,Aj27r- Heucc, wc get the following inequality: 

L'^£,Aj27r < ['^£+1/4 ,aJ 27r = ['^£+1/2 ,aJ 27r < ['^£+3/4 ,aJ 27r = L'^£+1,aJ 27r, 

which implies that a never goes below an integer multiple of 27r that it passes. This means 
ae > and ae — > ae for i > ii. 
Lemma [2U1 provides the bound 

|Ea£ — a^l < (a — a^) + c(y + v^) + Ed;f_i/2 + ^ 6jEdj. 

j=ei 

Since \E,ae — > Ea^, inequality ( l68ll follows for a = via Lemma [2T1 

In order to deal with the d = case in (l68|l we define T G Z/2 the first time when 

ttT > o^. Note that a can only pass an integer multiple of 27r in the i £ + 1/4 or 

£ + 1/2 — )'£ + Y4 steps, and (f evolves deterministically in these steps. This means that 

T — 1/2 is a stopping time with respect to J^j, j G Z/2. 

For large enough /C we have j^^^ < ^ • Then by Lemma [11] we get the uniform bound 

aj+i/4,A - ai,A < criQ , J < ^2, J e Z/2. 

By the strong Markov property and the bound flS^ we get 

E[(ar - a^)l(T < £2)] < cng and E[(ar+i/2 - a^)l(T < £2)] < ctiq 

Using this together with the first part of the proof and the strong Markov property again 
we get 

E{ai, - a^)+ = E [l(T < £2)^ [ae, - Jr]] 

< ci(E[(aT-a^)l(r<£2)] + + y) 

< cUV^ + y), (69) 



Lemma |20] gives 



i-2 

\Eai - a^\ < (a^ - a) + c(y + v^) + Ed^_i/2 + ^ 6jEdj. 
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Then by (169!) and the identity \a\ = —a + 2a"'' we get 

E\ae-a^\ < \Eae -a^\+ 2E(a^ - a^)+ 

e-2 

< (a^ - a) + c{y + v^) + Ed^-i/2 + ^ 6j-Edj. 
Since ag < — a^|, the Gronwall-type estimate in Lemma [2T] imphes fl68|) with a = a^. □ 
6.2 The uniform limit 

Proposition 22. Assume that K, = /C„ wi/i /C — oo anc? i/iai tiq "'^/C(nJ^^ VI) — t- 0. Then, 
(Pn2,o modulo 271 converges in distribution to the uniform distribution on (0,27r). 

Proof. We can use exactly the same argument as in Proposition 33, [13]. We show that 
given e > 0, every subsequence has a further subsequence along which fn2,o modulo 27r is 

1 /3 

eventually e-close to the uniform distribution. We set ^ = [/C(?^i V 1)J and pick r = r(e) 
with tC, < n2. Because of nQ^K,{n\^^ V 1) — )■ we will be able to let r — )■ oo. We will show 
that for any fixed r the distribution of v^nj.o ~V^n2-r5,o given Tn^-ri is asymptotically normal 
with a variance going to oo as r — )■ oo. From that the statement will follow. 

Note that the arguments of Proposition [16] can be repeated for the evolution of (y9„2_^g+^^o 
with < £ < which gives that v5n2,o conditioned on Tn^-T^ converges to (y9o(l — (1 + r)^-*^) 
where v9o(t) is the solution of (|5T]) with A = 0. This is just a normal random variable, 
its variance is given by the integrating the sum of the squares of the independent diffusion 
coefficients on [0, 1 — (1 + t)^^]. This is at least as big as the variance coming from the dW 
term which gives jj^ '^^^''^ I3{i~t) '^^^^ to oo if r — )■ oo as required. □ 

7 Last stretch 

The purpose of this section is to prove that on the interval [^2,?^] the relative target phase 
function af^ = ipf^ — ipf^ does not change much. 

Proposition 23. For any fixed A G M and IC > we have a®^ — ^ 0. 

1 /3 

The length of the interval [n2,n] is equal to ni + lC{n^ V 1), up to an error of order 1. 
By taking an appropriate subsequence of n (see Remark [H]) we may assume that Ui has a 
finite or infinite limit. We will consider these two cases separately. 
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Proof of Proposition in the limrii < oo case. By (fT2|) we have that |m — n — con- 
verges, by taking an appropriate subsequence we can assume that the hmit also exists without 
the absolute values. Note that condition f lTOj) implies that Ymvm/n > 1 and thus m — n — )■ oo. 

We may also assume that n — n2 is eventually equal to an integer S,- From there we 
proceed similarly as in [13]. We first note that by f l33|l and (|36|1 we have 

where R^,a = M^^J^^M^^J^^ From ([28]) and (|29]) we get 

T-, '^^n-j ~^ f^n-jYn-j ~ S„_jA 

r.R„_j- A - 



Pn—jPn—j+1 

where A = /i„ + A/(4ny^). Using (IT^ . = a/j — 1/2, = i^m — ri+j — 1/2 and 

P 

m — n — > oo we get that r.R„_j^A — ''^•Rn-i.o — ^ for any fixed j. This leads to 



a® ,A = ^ti,X*Qn-i-l - ^®_5,o*Q;-5-i ^0. □ 



If limni = oo then we will need the edge scaling results proved in [5] which are sum- 
marized in Theorem [3] of the introduction. The initial condition /(O) = 0, /'(O) = 1 for 
the operator "H^ in the theorem comes from the fact that the discrete eigenvalue equation 
for A„^m^n,m with an eigenvalue A = {^/n + y/m)^ + ^^^^^~06~~^ is equivalent to a three- 
term recursion for the vector entries w^^a (cf. f ll6p ) with the initial condition w^^^ = and 

By [9], Remark 3.8, the results of [9] extend to solutions of the same three-term recur- 
sion with more general initial conditions. We say that a value of u is an eigenvalue for a 
family of recursions parameterized by u if the corresponding recursion reaches in its last 
step. Suppose that for given ( G [— oo, oo] the initial condition for the three-term recursion 
equation satisfies 

{ro,u - 1) — ^ C, (70) 



uniformly in v with tq := Wi^u/wq^^. Here the factor ^^^^^^^^3 is the spatial scaling for the 
problem ([9], Section 5). Then the eigenvalues of this family of recursions converge to those 
of the stochastic Airy operator with initial condition /(0)//'(0) = C,. The corresponding 
point process is also a.s. simple and it will satisfy the following non-atomic property: for 
any a; G M we have P(x G S<^) = (see [9], Remark 3.8). Similar statement holds at the 
lower edge if liminf m/n > 1 with (r„ i, + 1)~^ in (ITH]) . (In this case one first multiplies the 
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ofF-diagonal entries of An^rn^m ^"7 ~1 before applying the arguments of [9], this will not 
change the eigenvalues.) 

If m/n — )■ 7 G [1, oo) then we can rewrite f lTOjl jointly for the upper and lower soft edge 

as 

^^^^^,n-^'\r,,Tir^C. (71) 

The multiplier is 1 in the case 7 = oo and it is always a finite nonzero value unless 7 = 1 
and we are at the lower edge. 

Proof of Proposition if lim ni = 00. By taking an appropriate subsequence, we may as- 



sume that /i„ — i/m — n is always positive or always negative. According to the proof of 
Lemma 34 in [13] we need to consider the family of recursions 



with initial condition 

for a given x G M and show that the probability of having an eigenvalue in [0, A] converges 
to as n — i- 00. 
We introduce 

n* = n — n2, m* = m — n2- 

Note that the recursion is determined by the bottom (2n*) x (2n*) submatrix of 

(!26|) which has the exact same distribution as the matrix A^l"^, . Thus we can consider the 

eigenvalue equation for A^!"^, with generalized initial condition 

(We will use * to denote that we are working with the smaller matrix.) We will transform 
this back to an eigenvalue equation for An*^m* and then An* ^m* A!^* m* with a generalized 
initial condition. 

Recall the recursion fl28l) and that r| = u}j^^/v}, f} = vl/u}. From this we get 



-J- = rlf: = rX ^ H IH ^ ' 

^0* ° ' '\ r*o p*o Pi) V Pi 



4 + mpn2)x + 5?(p„j)4 ) #^ (72) 



PO PO/ 
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where A = + Denote the solution of the generahzed eigenvalue equation for A^* 



,m* 



corresponding to A with ul,vl, .... If we remove the conjugation with D^"-*) from 
then from (172|) we get 

|^(-| + (a(p,.,). + »(p„J)A)M|. (73) 

By Remark [7] this is exactly the initial condition for the generalized eigenvalue equation for 
An*,m*Al.^^. with eigenvalue = /i^ + ^ + 

The spectrum of the matrix An*^rn*^*m* is concentrated asymptotically to [(v^m* — 
a/?2*)^, (Vm* + \fn*Y]. A direct computation shows that ^Jrrll ± = where we 
have + if /i„ — — n > and — otherwise. This means that if yU„, < -^/m — n then our 
original bulk scaling around /i„ corresponds to the lower edge scaling of An*^rn*A'^* and if 



fin > \pm^-n then we get the upper edge scaling. (Note that because of our assumptions if 



/i„ < v^m — n than liminf m/n > 1.) 

Again, by taking an appropriate subsequence, we may assume that m* /n* — )■ 7* G [1, C)o]. 
We first check that the initial condition ( I73|) satisfies ( I7T]) . this is equivalent to showing that 
=F 1^ converges in probability to a constant. Since ni — ?■ 00 we have 

m* >n* ^ 00, pn^ = ±1 + iv^(n*)-^/^ + o{{n*)-^'^), 
m = n* - lC{n*f''\ mi = m* - /C(r^*)^/^ 

Since 1 as £ — 00 we get that in probability 

hm (r2*)i/3 T 1) = v^x (1 ± (7*)-'/') 

and the convergence is uniform if we assume that A is bounded. This means that we may 
apply the edge scaling result, we just need to convert the scaling from the bulk to the edge. 
In the bulk we were interested in the interval I = [/x^, (/i„, + Atiq ^^^/4)^]. If we apply the edge 



scaling, then this interval becomes f-r!^ r^ TUsi^ ~ {^/m* ± a/ti*)^). Using our asymptotics 
for fin we get that for the end point of the interval 

ifii - iv^±v^r) ^ -^(1 ± (7*)-^/^)^/^ 



and for the length we get 



lim (fin + An. '^VA? - fil ' , = ^ ll ' lim 



1/6 (^*)l/6 X 



(v^±V^)^/3 (y^± 1)1/3 
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(We used that if 7* = 1 then we must have + in the ±.) This means that after the edge 
rescahng the interval shrinks to a point, meaning that the probabihty that our original 
recursion has an eigenvalue in [0, A] converges to the probability that the limiting edge point 
process has a point at a given value which is equal to 0. □ 



8 Appendix 

In this section we provide the needed oscillation estimates. 

Lemma 24. Suppose that 2tx > 61 > 62 > ■ ■ ■ > 0^ > ^ and let si = Yl^j=i ^j- Then 

m 



Proof. The first inequality is the same as Lemma 36 from [13] and the second inequality is 
straightforward. □ 



Lemma 25. Let 1 < ii < £2 ^ no and F^^^^^ = Ylf=£i Vj fof i = 1,2 and g£ ^ C . Then for 
i = 1,2." 



£2-1 



£2 



£2-1 

^ iT^iuiildiiW + 5Z \^eii/2\\9e+i - 9e\ 



We also have the following estimates: 

1^(1) I 



< c{l + n/^k2^^^) 



(74) 
(75) 



(76) 



|ir(2) I < 



{\pIpI + + l + n'/'k-'^') if (^) or (^^) 
IpI + + 1 + n'/^k^'^^] if f^**) 



where the conditions are given by 



'm 


— n 


m 


— n 



[-k-k-k) : /i„ < \/m — n with /cq, k\ < y/mifii. 



(77) 
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Proof. The bounds (17^ and (175|) follow from partial summation. In order to prove the 
bounds on F^^\^ we will apply Lemma |2U but we need to consider various cases. Note that 
the constant c might change from line to line. 



Case 1: /i„ > ^/■m — n, 3?p^ > 0. 
We have the bounds 

k^/^{ni + ky^l"^ < argp^, n^'^{ni + A;)"^/^ < 7r/2 - argp^, 
k^^'^{mi + k)~^^'^ < argp^, m|^^^(mi + k)'^^"^ < tx 12 — argp^. 

and a.Tg{piPi) is decreasing. The sequence a.Tg{pjpj) satisfies the conditions of Lemma 
and using mi > rii and mi > cn > cuq we get the bound 

\J2vi\<c ((m + k2Y/''k-'^^ + {ki + mi)i/2^-'/') < cin]'^-^'^ + 1). 

We have 

y/miTii — k i\/k y/n^ 
^^^^ \/k + m\\/k + ni y/k + mi^/k + rii 

which means that if fci , /c2 > -v/mirii or A;i , A;2 < ^mini then we can use Lemma [2l] for 

the sequence a.Yg{pjpf). (In the first case 7r/2 < arg(p^p^) < vr while in the second case 
< aig{piPi) < 7r/2.) From the lemma we get 

and explicit computation together with mi > cn > ck gives 



iVfc [^/rnl + y/ni) 

This finishes the proof of the lemma in this case. 
Case 2: /i„ < y/m — n, ^pe < 0. 
Now we have 

k + ^/m^y/n^ iVk {y/m[ — y/n^) 



Plk = rr^ = H rr^ ^- C^^) 

VA; + rraiV«^ + '^i \/k + miy/k + ni 

meaning ir /2 < axg{pipi) < it. Differentiation of the real part shows that axg pipi decreases if 
k > y/miUi and then it increases. This means that we can apply Lemma [241 for the sequences 
a.Tg{pjpj) or a.Tg{pjpj) if /ci,/c2 > y/nimi and the reversed versions of these sequences if 
/ci, /c2 < y/n^mi. From (I79|l we get 



TT < argp2p2^ 27r - argp2p2 > 2^^^^^^1 > 0(^^1/2^^/' + 1)^1 (80) 
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where we used — n\ > c^frn{ which follows from mi > cn > crii. Applying Lemma |2 
to aig{f>1pf) (or the reversed sequence) we get 



c(n}^^/c^ ^^"^ + 1) if ki,k2 > y/riirrii, 



t=ti 



c(n 



\^'^k2 + 1) if ki,k2 < ^n\m\. 



Noting that k^^^"^ < k2^^'^ this proves fITH]) in Case 2 ii ki,k2 < ^Jn^m{ or k\,k2 > ^n\m\. 
If k\ > y/riinii and k2 < y/riimi then we can cut the sum in to parts at y/n^m^ and since 
for k > y/mirii we have n\^'^k~^^'^ < 1 we have (1761) in this case as well. 
To prove (1771) we apply Lemma [241 to aTg{pfpf) (or its reversed) to get 



IE"? 

£=£i 

From this (I77p follows by noting that 



c {\pIpI + 1| ^ + IpIpI - 1| ^) if ki, k2 > y/n^, 

c {ipIpI + ir^ + ipyi - ir^) if ki, k2 < y^^. 



2V k [^/ml — ^/ni) 

where the first equality is explicit computation and the inequality is from ( IHOj) . □ 
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